UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


MECHANICAL  BEHAVIOR  AND 
FATIGUE  STUDIES  OF  RUBBER 


COMPONENTS  IN  ARMY 
TRACKED  VEHICLES 

H.R.  Brown,  J.L.  Bouvard,  D.  Oglesby,  E.  Marin,  D.  Francis,  A. 
Antonyraj,  H.  Toghiani,  P.  Wang,  M.F.  Horstemeyer,  M.P.  Castanier 


This  project  was  funded  by  SimBRS  and  performed  at  the  Center  for  Advanced 
13  August  2010  Vehicular  Systems  (CAVS)  at  MSU 


GVSETS 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

13  AUG  2010  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Mechanical  Behavior  and  Fatigue  Studies  of  Rubber  Components  in 
Army  Tracked  Vehicles 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

W56  HZV-08-C-0236  (SimBRS) 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

H.R.  Brown;  J.L.  Bouvard;  D.  Oglesby;  E.  Marin;  D.  Francis;  A. 
Antonyraj;  H.  Toghiani;  P.  Wang;  M.F.  Horstemeyer;  M.P.  Castanier 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

The  Center  for  Advanced  Vehicular  Systems  (CAVS)  at  Michigan  State 
University 

8.  PERFORMING  ORGANIZATION  REPORT 

NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

US  Army  RDECOM-TARDEC  6501  E  11  Mile  Rd  Warren,  MI 
48397-5000,  USA 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

TACOM/TARDEC 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

21108RC 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Presented  at  NDIAs  Ground  Vehicle  Systems  Engineering  and  Technology  Symposium  (GVSETS),  17  22 
August  2009, Troy,  Michigan,  USA,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION 

OF  APSTR  apt 

18.  NUMBER  19a.  NAME  OF 

HF  PAfTRQ  PFQPnMQTFT  P  PFDQflM 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  S  AR 

unclassified  unclassified  unclassified 

39 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


•  Monotonic  approach 

-  Internal  state  variable  (ISV)  model  for  metals 

-  ISV  modeling  strategy  moved  to  glassy  polymers  (Bouvard  et  al.,  2010) 

-  Current  efforts  to  apply  ISV  modeling  strategy  to  elastomers 

•  Fatigue  approach 

-  Researchers  have  historically  separated  fatigue  crack  initiation  and  propagation 

-  (McDowell  et  al.,  2003)  refined  the  earlier  crack  stages  into  incubation  and 
microstructurally  and  physically  small  crack  growth,  greatly  increasing  accuracy 

-  Microstructure  has  been  incorporated  into  the  multistage  modeling  for  metals  at 
CAVS 

-  Researchers  have  typically  only  investigated  long  crack  for  elastomers  (Mars  and 
Fatemi,  2003;  Busfield  et  al.,  2002;  Chou  et  al.,  2007) 

-  Current  efforts  are  to  add  MSC/PSC,  INC  to  fatigue  modeling  of  elastomers  and 
incorporate  microstructure 
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Overview 


Experiments  and  Stereology  to  Capture 
Structure-Property  Relationships 


— 0.1/S 

— 0.01/s 
— 0.001  Is 


System  Level  Modeling 
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Multi-Scale  Internal  State  Variable 
Modeling 


Material  Database  and 
FEA  Material  Model 
Development  for  Polymer 
Materials 


Component Level 


Validation 


Material  Level 
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•  Components  of  focus  on  tank  track 

-  Bushings 

-  Road  wheels 

-  Suspension 

-  Backer  pads 

•  Extreme  Loading  conditions 

-  High  temperature 

-  High  friction 

-  Complex  loading 

•  Road  wheel  backer  pad  failure  at 
one-half  of  the  design  target  mileage 
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Macroscale  MSU  ISV/MSF 
Models  Implementation  and  Use 


Multiscale 

Materials 

Modeling 


Physics  Validation 
And 

Numerical  Verification 


t 


boundary  conditions 
loads 

temperature 
strain  rate 
history 
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ISV=lnternal  State  Variable 
MSF=MultiSteg&Fatigue 


Note:  model  can 
be  implemented 
in  other  FE  codes 
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“  J  Multiscale  Experiments 


I 


1 .  Exploratory  exps 

2.  Model  correlation  exps 

3.  Model  validation  exps^ 


Nanoscale 


Model 

Cohesive  Energy 
^Critical  Stress 


Analysis 

Fracture 
Interface  Debonding 


FEM  Analysis 

Idealized  Geometry 
Realistic  Geometry 
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Material:  Styrene  Butadiene 
Rubber 


H 

V 
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•  Random  copolymer  -  styrene  and  butadiene  are 
randomly  distributed  throughout  the  polymer  chain 

•  3:1  butadiene  to  styrene  by  weight 

•  Commercially  used  in  a  wide  range  of  projects 


General  properties  of  SBR 


Property 

Value 

Units 

Glass  Transition  temperature 

-40 

°  c 

Temperature  range 

-28  to  76 

°  c 

Tensile  strength 

4.8 

MPa 

Stretch  Limit 

150 

% 

Density 

91.5 

Lbs/c u  ft 
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styrene  unit  butadiene  unit 

©1997  Encyclopaedia  Britannica,  Inc. 
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Experimental  Methods:  DMA 


Dynamic  mechanical  analysis 

-  TA  Instrument  Q900  Dynamic  Mechanical  Analyzer 

-  Rectangular  1 .5  in  x  0.1 35  in  x  0.06  in  specimen  with  DMA  tensile 
clamps  (ASTM  D4065-01 ) 

-  Oscillated  at  1  Hz  for  a  range  of  temperatures  to  include 
temperature  transitions 

-  Tg  measured  using  midpoint  method 


13  August  2010 


GVSETS 


Experimental  Methods:  Set-up 


MSTV 

MODELING  rho  SIMULATION,  TESTING  rkd  VRLIDHTION 


w 


Compression  Test 


Hopkinson  Bar 


Tension  Test 


Low  strain  rates 
(1  O'4  S'1  to  10-1  s-1) 


Low  strain  rates 
(1  O'4  S'1  to  10-1  s-1) 


Specimen  Geometry 


0.5  in 


0.25  in 


Geometry  of  0.5  L/D  ratio 


Low  frequency 
(2Hz) 


Scaled  from  ASTM  D41 2 
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Experimental  Methods: 
Monotonic  Loading 


Stress  state  dependence 

-  Tension 

•  Strain  controlled  from  extension  to  get  local  strain  rate 

•  Strain  in  the  gage  measured  by  laser  extensometer 

-  Compression 

•  Extensometer  mounted  on  platens  to  remove  compliance 

•  Strain  controlled  from  extensometer 

Rate  dependence 

-  Strain  rates  of  0.001 , 0.01 ,0.1,  1200,  2000,  2600,  3000  /s 

Temperature  dependence 

-  -5C,  23C  and  50C 
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Experimental  Methods:  Fatig 
Loading 


•  Uniaxial  tension 

•  Servo-hydraulic  load  frame 

•  R  =  0.5,  freq  =  2  Hz 

•  A£/2  =  20.3,  22.3,  28.1,31.9 
and  36.3% 

•  Displacement  control 

•  Displacement-strain 
correlation  made  using  laser 
exten  so  meter 
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3.5 


Grip  to  Grip  Diplacement  (mm) 
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Results:  DMA 


•  Storage  modulus  curve 
shows  two  transition 
temperatures 

•  First  corresponds  to  the  glass 
transition  temperature  (Tg)  at 
-40°  C 


•  Second  at  -10°  C  possibly 
due  to  fillers 


Macroscale 

Experiment 

Uniaxial/torsion 
Notch  Tensile 
atigue  Crack  Gro 
Cyclic  Plastuciiy^^ 


ntinuum 


Cyclic  Plasticity 
Damage. 

hK  M  Analysis 

Torsion/Comp 

Tension 
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Tan  Delta 


Results:  Rate  Dependence 
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0.1  Is 
0.01 Is 


2  0.4  0.6  0.8  1 

True  Strain  (mm/mm) 


True  Strain  (mm/mm) 


50 


True  Strain  (mm/mm) 


Rate  dependence  -  material  stiffened  with  increasing 
strain  rate 


Compression  -  strong  hysteresis  and  residual  plastic 


L  4 


j 


Results:  Temperature 


j 


Dependence 
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9 


True  Strain  (mm/mm) 


50 


True  Strain  (mm/mm) 


•  Temperature  dependence  -  material  softened  and  had 
an  increase  in  strain  to  failure  at  higher  temperatures 
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^  Results:  Microstructure  for 
Monotonic  Loading 
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•  Particle  debonding 

-  Major  mechanism  of  failure 

-  0.5  -  200  |jm  particles  debonded  during  deformation 

-  Large  particles  are  (1)  agglomerations  of  aluminosilicate  (clay) 
and  (2)  ground  calcium  carbonate 
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MSU-EMC 


SEI  5.0kV  X900  WD  9.9mm  10/um 
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Results:  Microstructure  for 
Monotonic  Loading  - 1  =  0.1/s 


» 


MODELING  rho  SIMULATION,  TESTING  rno  VALIDATION 


'  '  '  ‘i "P 

■ 

WL .. 
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'  ,  *v  *  •  .  .  * 


Room  temperature  tension  at  t=  0.1 /s 
Specimen  A  failed  at  a  lower  strain  than  specimen  B 

Specimen  A  showed  a  slightly  weaker  stress-strain  response  than  specimen  B 

Due  to  the  large  agglomerates  of  undispersed  aluminosilicate  and  particles  debonding  from  the  matrix 
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Results:  Microstructure  for 
Monotonic  Loading  - 1  =  0.01/s 


If 

V 
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•  Room  temperature  tension  at  i=  0.01  /s 

•  Specimen  A  failed  at  a  lower  strain  than  specimen  B 

•  Due  to  the  large  agglomerates  of  undispersed  aluminosilicate  and  particle  debonding 
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Results:  Microstructure  for 
Monotonic  Loading  - 1  =  0.001/s 


If 

V 
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•  Room  temperature  tension  at  i=  0.001 /s 

•  Specimen  A  failed  at  a  lower  strain  than  specimen  B 

•  Due  to  the  large  agglomerates  of  undispersed  aluminosilicate  and  particle  debonding 
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Macroscale  MSU  ISV/MSF 
Models  Implementation  and  Use 


Microstructure  and 
experiments 


Multiscale 

Materials 

Modeling 


Physics  Validation 
And 

Numerical  Verification 


boundary  conditions 
loads 

temperature 
strain  rate 
history 
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ISV=lnternal  State  Variable 
MSF=MultiStage  Fatigue 


Note:  model  can 
be  implemented 
in  other  FE  codes 
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Large  variety  of  models  exists  for  polymers:  Krempl 
(1995),  Tervoort  (1998),  Boyce  et  al.  (1988),  Richeton  et 
al.  (2007),  L.  Anand  et  al.  (2009),... 


Rheological 

representation  of  the  constitutive  model 


/y //////// 


Model  generally  used  for  polymers  are:, 
>  based  on  spring/dashpot 


Linear  spring 

Nonlinear  viscous  i  _j_  I 
dashpot  ■— | — * 


L  _ 


Network  (A) 


Hyperelastic 
rubbery  spring 


Network  (B) 


o  = 


+ 


B 


Hierarchical  Multiscale  Approach 

Development  of  ISV  material  model: 

-  Kinematics 

-  Thermodynamics  select  physically-based  ISVs 
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ISV  Model:  Extension  to 
Elastomers 


Jr 
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Regime  I:  Hyperelastic  mechanism  induced  by  bond  stretching/rotation 
Regime  II:  Strain  hardening  induced  by  crosslinking,  entanglements,  and  particles 


Regime  III:  Chain  alignment  and  chain  stretching  between  crosslinking  and  possible  chain  crystallization 


MchipnChi 
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Amorphous  Polymers 


■  ■  K£J  m  w 
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F  =  FeFp 


Tj^e  :  elastic  mechanisms  such  as  bond  stretching  and  chains  rotation/torsion 
inducing  the  different  conformations  of  the  intramolecular  structure 

Fp  :  time-dependent  inelastic  mechanisms  such  as  permanent  chains  stretching  and 

rotating  but  also  the  dissipative  mechanism  due  to  the  relative  slippage  of  molecular 
chains 
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ISV  Model:  Kinematics  and 
Thermodynamics 
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Multiplicative  decomposition  of  deformation  gradient 
F  =  FeFp  Deformation  gradient 

Kroner-Lee  decomposition 

r  =  det  Fe ;  Jp  =detFp; 

Kinematics  ->  velocity  gradient  1: 

1  =  FF  1 


1  =  le  +  FeLpFel 

•  le  =  FeFe  1 ;  le  =  de  +  we 

.  Lp  =  FPFP~‘  LP=DP  +  WP 


Decomposition  into  symmetric  and  skew  part 


•  Assumption  on  plastic  flow 

•  Jp  =  1  Flow  is  incompressible 

.  wp=o  Flow  is  irrotational 

•  Thermodynamics 

Clausius-Duhem  inequality  (Gurtin  and  Coleman^  1967;  Gurtin  and  Anand,  2003) 

D,  .  r  ^  S:FeT 

—  jje  \|/d  v  <  J  <r :  1  d  v  local”Tegrals  \ 

L>t  B _ J3  S_=  J 
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d  F  +  Mj  Dp  —  \j/  >  0 

M  =  CeS 
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e-1 


e 1aFe T 
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>■'  ISV  Model:  Kinematics  and 
Thermodynamics 
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Helmholtz  Free  Energy: 

¥ 


:  Strain  field  induced  by  internal  strain  field  related  to  the  intermolecular  chain 
interaction  in  Regime  I  and  II  (van  der  Waals  forces  mainly) 

^2:  Strain  field  induced  internal  strain  field  related  to  the  crystallization  in  Regime  III 

:  Stretch-like  tensor  giving  direction-dependent  (kinematic)  hardening  effect  induced 
by  the  chain  stretching  between  entanglements/crosslinking  at  large  strain 

+ 

dCe  1  d^2  2  ap 


Clausius-Duhem  Inequality: 


f 


S-2-^ 


da 


FeTdeFe  +  M  :  Dp 


A 


.X+d~i-x+8-i-.t 


St, 


K1 


dt>  2 


K. 


ap 

a 


>0 
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ISV  Model:  Summary 
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Elasticity:  Cauchy  Stress 

F  =  FeFp ;  Fe  =  ReUe ;  Fe*  =(je)"‘/3Fe;  be*=Fe*F 
T  =  FeSFeT;  T  =  [1B  (lev(bc* )+  K B (je)l ;  o  =  J^t 
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Parameters  {m  ^p  >  y,  ti0,  <^0,  ^sat,Rsl,RS2,cKi  ,ca,x,L,[xR  J 
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ISV  Model:  Model  to  Experii 
Comparsion 


ISV  model  predicts  loading  path, 
unloading  path  and  time  dependence 
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0  0.05  0.1  0.15  0.2 

True  Strain  (mm/mm) 


ISV  model  response  under  cyclic 
loading  at  1  Hz.  It  shows  significant 
hysteresis  in  the  first  cycle  and  cyclic 
relaxation 

26 
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Microstructure  and 
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Multiscale 

Materials 

Modeling 


Physics  Validation 
And 

Numerical  Verification 


boundary  conditions 
loads 

temperature 
strain  rate 
history 


13  August  2010 


ISV=lnternal  State  Variable 
MSF=MultiSteg&Fatigue 


Note:  model  can 
be  implemented 
in  other  FE  codes 
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Results:  Fatigue  Life 


•  Power  law  fit:  0.6587^  )~0093 
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Results:  Fatigue  Loading 


If 

V 
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•  Significant  cyclic  stress  softening  occurred  for  all  strain 
amplitudes 
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True  Stress  (MPa)  True  Stress  (MPa) 


Results:  Hysteresis 
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Ae/2  =  0.36 


Ae/2  =  0.28 
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A£/2  =  0.32 


A£/2  =  0.22 
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Results:  Hysteresis 


0  0.2  0.4  0.6  0.8 

True  Strain  (mm/mm) 


1.2 


Half  cycle  at  strain  amplitudes  of 
0.22,  0.28,  0.32  and  0.36. 
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to 

to 

0) 


CO 

Q) 

D 


Last  cycle  at  strain  amplitudes  of 
0.22,  0.28,  0.32  and  0.36. 
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Results:  Microstructure  for 
Fatigue  Loading 
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Particle  debonding 

-  0.5  -  200  |jm  particles  debonded  during  deformation 

-  1 00  -  200  pm  particles  initiated  significant  fatigue  cracks 

-  2  particles  in  length  scale  of  focus:  calcium  carbonate  and 


agglomerations  of  aluminosilicate  (clay) 
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Results:  Microstructure  for 
Fatigue  Loading 


If 

V 
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Results:  Microstructure  for 
Fatigue  Loading 


If 

V 
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0.5 


0.3 


0.2 


0.1  - 


icr  ioJ  io4  io&  ioD 

Cycles  to  Failure,  Nf 
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&  Results:  Microstructure  for 
Fatigue  Loading 
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Multistage  Fatigue  Model 


Mi 


'JFf 
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NtotarNinc+ N  MSC+ ^  PSC+ N  |_C 


Incubation 


Ay. 


max 


=  Cinc N“c  Cmc  =  Cn  +  z(Cm  -  C. )  and  C.  =  0.24(l  -  (*)) 


MSC/PSC  Growth 


^  da  ^ 


\dN  j  MSC !  psc 


ACTD  =  f((/>)Cn  (l) 


=  G(ACTD-ACTDth ) 


Initial  crack  size 
a,  =  0.5625 D, 


UAa 


a 


+  C/(1) 


A  y 


P  \ 
max 


T 


J _ I  L 


/ 


T 


macro  j 


/(^)  =  1  +  co\  1  -  exp 


HCF  loading  dominated  LCF  loading  dominated 

f  i  A' 


Ad  =  2# 


,  -|0.5 

3  Act.  Affj 


2  2  2 

^  da  ^ 


Porosity  term 

V  ^th  )  J 

+  (l-^)A<r1  Multiaxial  term 


ACTDth  =nonpropagating  crack  threshold 
1 


t/  = 


1-7? 


Mean  stress  term 


LC  Growth 


\dN  j  LC 
r  da  ^ 

\dN  JLC 


=  a{t-t0)  for  T0<T  <ta 


at  >  250  pm,  or 


=  BTfl 


for 


ta<t<tc 


da 

dN 


=  max 


da 

\dN  j  Msc/psc 


da  \ 


dN 


J  LC 
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•  DMA  testing  was  performed  to  investigate  the 
viscoelastic  properties  and  transition  temperatures 

•  Material  exhibited  time  and  temperature  dependence 

•  Debonding  of  calcium  carbonate  particles  and 
aluminosilicate  agglomerates  on  the  order  of  50  to  200 
pm  lead  to  specimen  failure  for  monotonic  loading  and 
initiated  fatigue  cracks  under  fatigue  loading 

•  The  ISV  model  captures  both  loading  and  unloading  as 
well  as  rate  dependence 
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•  The  MSF  model  equations  need  to  be  extended  to 
elastomeric  materials  as  well  as  calibrated  and  validated 
on  SBR 

•  Fatigue  experiments  need  to  be  conducted  at  lower 
strain  amplitudes  to  investigate  the  high  cycle  fatigue 
response  of  the  material 
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